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Abstract 

The barriers standing against the formation of superheavy elements and their con- 
secutive a decay have been determined in the quasimolecular shape path within a 
Generalized Liquid Drop Model including the proximity effects between nucleons 
in a neck, the mass and charge asymmetry, a precise nuclear radius and the shell 
effects given by the Droplet Model. For moderately asymmetric reactions double- 
hump potential barriers stand and fast fission of compact shapes in the outer well is 
possible. Very asymmetric reactions lead to one hump barriers which can be passed 
only with a high energy relatively to the superheavy element energy. Then, only 
the emission of several neutrons or an a particle can allow to reach an eventual 
ground state. For almost symmetric heavy-ion reactions, there is no more external 
well and the inner barrier is higher than the outer one. Predictions for partial a 
decay half-lives are given. 
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1 Introduction 



Using heavy-ion reactions of mean asymmetry (Ni on Bi and Zn on Pb 
[1,2]) and more recently highly asymmetric reactions (Ca on U, Pu and Cm 
[3,4,5]) new very heavy elements have been synthetized. In these very long 
and difficult experiments the number of detected events (a decay chains) is 
very low and the interpretation of the primary experimental data is discussed 
[6,7,8]. Actually, it is assumed that the lower limit for the fission barrier heights 
of these heaviest elements is around 6 MeV [9]. To form other superheavy 
nuclei, many other reactions covering all the asymmetry range are planned or 
are proposed using other isotopes, inverse kinematics and radioactive nuclear 
beams (and targets) [10,11,12]. 

The purpose of the present work is, firstly, to study the potential barriers 
standing in the quasimolecular shape path for these new suggested entrance 
channels and to look at their dependence on the next proton magic number 
(114, 120 or 126). The potential barriers relative to already realized experi- 
ments have been determined in a precedent paper [13]. The second objective is 
to determine the partial a decay half-lives of these superheavy elements from 
the Q values given by a recent version of the Thomas- Fermi model [14]. 

The Generalized Liquiq Drop Model which has allowed to describe most 
of the fusion [15] , fission [16,17], light nucleus [18] and a emission [19] data 
has been used once again. The shell effects have been derived from the Droplet 
Model formulae [20]. 



2 Quasimolecular shapes 



To describe the entrance channel leading from two spheres to one sphere 
and the alpha emission a shape sequence leading rapidly to the formation of a 
deep neck while keeping almost spherical ends has been retained (see Fig. 1). 
Mathematically, these quasimolecular shapes correspond to two joined elliptic 
lemniscatoids [15] assuming volume conservation. Analytic formulae are avail- 
able for the main shape- dependent functions : volume, surface, distance be- 
tween mass centres and moment of inertia. Similar compact shapes have been 
also adopted to determine the energy of superheavy nuclei [21]. For the exit 
channel these shapes implicitly suppose that the a particle is rather formed 
at the surface of the parent nucleus. 
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3 Generalized Liquid Drop Model 



The macroscopic energy is determined within the Generalized Liquid 
Drop Model previously defined and used to describe the fusion and fission 
processes as well as the light nucleus and a emissions [15,16,17,18,19]. 

E = E Vo l + E Su rf + E Coulomb + E Prox . (1) 



When the colliding nuclei or the spherical fragments are separated: 

E Vol = -15.494 [(1 - l.d>ll)A l + (1 - 1.8/ 2 2 )A 2 ] MeV, (2) 

Esurf = 17.9439 [(1 - 2.QI 2 )A 2 / 3 + (1 - 2.QI%)A 2 2 /3 ] MeV, (3) 
Ecouiomb = O^Zl/R, + 0.6e 2 Z 2 2 /R 2 + e 2 Z 1 Z 2 /r, (4) 

where A iy Z iy Ri and Li are the masses, charges, radii and relative neutron 
excesses of the two nuclei, r is the distance between the mass centres. 
The radii Ri are given by: 

Ri = (1.28A 1 / 3 - 0.76 + 0.8A; 1/3 ) fm. (5) 



This formula allows to follow the experimentally observed increase of the ratio 
ri = Ri/A 1 / 3 with the mass; for example, r = 1.13/m for 48 Ca and r = 
1.18/m for 248 Cm. 

For one-body shapes, the volume, surface and Coulomb energies are defined 

as 

E Vol = -15.494(1 - 1.8I 2 )A MeV, (6) 
E Surf = 17.9439(1 - 2.QI 2 )A 2/3 (S/ArtR 2 ) MeV, (7) 

Ecouiomb = O.Qe 2 (Z 2 /R ) x 0.5 J (V(6)/V )(R(6)/R ) 3 sin 9d6. (8) 

/ is the relative neutron excess and S the surface of the one-body deformed 
nucleus. V(6) is the electrostatic potential at the surface and V the surface 
potential of the sphere. The radius of the compound nuclear system has been 
calculated from the radii of the two fragments assuming volume conservation. 

The surface energy results from the effects of the surface tension forces 
in a half space. When there are nucleons in regard in a neck or a gap between 
colliding nuclei or separated fragments an additional term called proximity 
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energy must be added to take into account the effects of the nuclear forces be- 
tween the close surfaces in regard. This term is essential to describe smoothly 
the two-body to one-body transition and to obtain reasonable fusion barrier 
heights. It moves the barrier top to an external position and strongly decreases 
the pure Coulomb barrier. 

E Prox (r)=2 1 J <5>[D(r,h)/b]2nhdh, (9) 

where h is the distance varying from the neck radius or zero to the height of the 
neck border. D is the distance between the surfaces in regard and b = 0.99/m 
the surface width. $ is the proximity function of Feldmeier [22]. The surface 
parameter 7 is the geometric mean between the surface parameters of the two 
nuclei or fragments: 

7 = 0.9517^/(1 - 2.6/ 1 2 )(l - 2.6/|) MeVfm- 2 . (10) 

In this GLDM the surface diffuseness is not taken into account and the prox- 
imity energy vanishes when there is no neck or gap. The combination of these 
GLDM and shape sequence has allowed to reproduce the fusion barrier heights 
and radii, the fission and the a and cluster radioactivity data. 



4 Potential barriers governing the entrance channel 

The macroscopic potential barriers (deformation energy) standing in 
front of some heavy- ion combinations leading to 270 110 and 302 120 compound 
systems are displayed in Fig. 2. Some reactions using very exotic nuclei are 
not still feasible. Within this GLDM the nuclear proximity energy introduces 
an inflection in the curve and the outer barrier top corresponds always to an 
unstable equilibrium between the repulsive Coulomb forces and the attractive 
proximity forces for two separated nuclei. For the highest asymmetries the 
external fusion barrier energy is higher than the spherical system energy and 
a potential pocket appears after crossing the barrier. In this description, the 
nuclear system can reach a quasi-spherical shape with little excitation energy 
via a tunneling process. This does not prove the stability of the formed system. 
For the exit channel all the symmetric and asymmetric fission and a emission 
barriers must be investigated. With decreasing asymmetry the outer well pro- 
gressively disappears as the outer peak. It will be very difficult to reach the 
eventual ground state of the superheavy element via these almost symmetric 
reactions. Predictions from such barriers for four asymmetric and almost sym- 
metric reactions leading to 246 Fm had been compared in a previous paper [15] 
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with experimental data on evaporation residues and fast fission events and the 
agreement was very good. 



In ordinary fusion studies, it is only the preceding barriers which are 
taken into account. The macroscopic energy does not reproduce very accu- 
rately the Q value since the nonuniform distribution of single-particle levels 
leading to the shell correction term and pairing energy is not taken into ac- 
count as also Wigner energy... In Fig. 3, the difference between the theoretical 
Q value given by the Thomas-Fermi model [14] and the GLDM one has been 
empirically added at the macroscopic energy of the compound nucleus with 
a linear (in r) attenuation factor vanishing at the contact point. Such a de- 
scription is valid if the nuclear system has enough time to relax and built its 
shells and pairs during its descent to a quasi-spherical shape. Then the barrier 
against reseparation is higher and the external minimum disappears for the 
highest asymmetries. 

In the following the shell corrections proposed in the Droplet Model 
[20,16] will be used and the formulae are recalled here. 

E She a = E s s tTa-29 2 )e^. (11) 



The shell corrections for a spherical nucleus are given by 

E s s f*J e = 5.8 [(F(N) + F(Z))/(0.5A) 2 / 3 - 0.26A 1 / 3 ] MeV, (12) 

where, for Mj_i < X < Mi, 

F(X) = Qi (X - Mi_i) - 0.6(X 5 / 3 - M?H). (13) 

Mi are the magic numbers and 

q t = 0.6(Mf /3 - M?Ii)/(Mi - M^). (14) 



Within this algebraic approach it is possible to study the influence of the se- 
lected value for the highest proton magic number. For the two highest neutron 
magic numbers the values 126 and 184 have been retained. 

9 2 = (5R) 2 /a 2 . (15) 



The distortion 6a is the root mean square of the deviation of the nuclear 
surface from the sphere, a quantity which incorporates indiscriminately all 
types of deformation. The range a has been chosen to be 0.32r . 
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Potential barriers for three very asymmetric reactions ('warm' reactions) 
A8 Ca+ 24A Pu, 48 Ca+ 248 Cm and 50 Ti+ 248 Cm are shown in Fig. 4 assuming dif- 
ferent hypotheses. Due to the high asymmetry and, consequently, to a low 
Coulomb repulsion and proximity energy the barrier against reseparation is 
wide and high. There is no double-hump barriers. Even for a subbarrier tunnel- 
ing of 6 or 7 MeV in the entrance channel and even if the shells and pairs have 
not enough time to develop (dashed curve) before reaching the quasispherical 
shape the nuclear system has enough energy to reach this quasispherical com- 
pound system. The excitation energy relatively to the ground state energy is 
more than 30 MeV, allowing the emission of several neutrons or an a parti- 
cle. The different hypotheses on the proton magic number do not change the 
global predictions in the entrance path. 

Potential barriers governing the medium asymmetry reactions ('cold' re- 
actions) m Ni+ 208 Pb, 70 Zn+ 208 Pb, 76 Ge+ 208 Pb and 82 Se+ 208 Pb are displayed 
in Fig. 5. A wide macroscopic potential pocket due mainly to higher prox- 
imity energy and Coulomb repulsion (relatively to warm reactions) appears 
at large deformations. Whatever the microscopic correction assumptions are 
double-hump barriers begin to appear for the Ni on Pb and Zn on Pb reac- 
tions. Experimentally the beam energy was respectively 239 and 257 MeV. 
Consequently, these successfull experiments correspond to subbarrier fusion 
processes in our approach and the quasispherical system can be reached by 
tunneling even if the shells and pairs are not completely built. For the two 
last reactions the inner macroscopic hill is higher than the external barrier. 
Incomplete fusion and fast fission events in the external pocket are the main 
exit channels since the neck between the two nuclei is formed and exchanges 
of nucleons may occur. If the reorganization of the single-particle levels is very 
rapid then the value of the proton magic number begins to play some role. 
So an open question is whether at large deformations the nucleon shells can 
take form to stabilize the nuclear system before investigating a peculiar exit 
channel. The pre or post equilibrium nature of the evaporation process of the 
excess neutron is also crucial. 

The same comments are valid for the two cold reactions shown in Fig. 
6 leading to nuclei of charge 118 and 119. A value of 120 for the next proton 
magic number naturally would slightly increase the possibility to form such 
superheavy elements. For still heavier mass and charge and more symmetric 
heavy-ion reactions (see Fig. 7) the potential pocket between the two peaks 
progressively disappears. A tunneling effect in the outer peak does not allow to 
reach a quasispherical state. Only a subbarrier fusion through the inner peak 
could eventually lead to a compound nucleus if the microscopic contributions 
grow rapidly at the beginning of the fusion process and if the charge of the 
compound system is the next proton magic number. 

The symmetric decay barriers for 283 112, 292 116, 295 120 and 311 126 are 
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given in Fig. 8. The microscopic contributions play the main role. The height 
and width of the one hump fission barrier is governed by the proximity of a 
proton or neutron magic number. The barrier height still reaches several MeV 
and can reach 10 MeV if the neutron and proton shell effects add each other. 

The characteristics of the macroscopic potential barriers for the new sug- 
gested reactions [10,11,12] leading possibly to superheavy elements are shown 
in table 1. The Q value is extracted from [14]. The energy of the external peak 
relatively to the sphere energy is positive for Z\Z 2 lower than around 2500. 
The position of the external barrier corresponds to the contact point of the 
colliding nuclei when Z\Z 2 is higher than around 2600. There is no more real 
external macroscopic barrier for Z\Z 2 higher than around 3200. 

In the hope of reaching the most efficient beam energy, formulas giving 
the 1-dependent fusion barrier height Ef us and radius Rf us are proposed. 



E fus>l (MeV) = E fus>l=Q (16) 

1(1 + 1) 



+ 



02081 (A 5/3 + A 5/3 ) + a050MlA2 ( A i /3+ ^ /3 ) 2 ( 1 - 908+ ii%- a0857fa ( ZlZ2 )) 2 ' 



with 



2.1388^^2 + 59.427(^ /3 + A^ 3 ) - 27.07 ln( 1 f 3 lZ \ /z ) 



{A^ + A 2 ^)(2.97 - 0A2\n(Z 1 Z 2 )) 



E fus , l=0 (MeV) = -19.38 + — ^ 1/3 „ — - 1 2 , (17) 



RfusAfm) = (A\ /3 + A L 2 /0 ) [1.532 + 6.48/^ 



l/3x 



(18) 



- o.ooo507(a} /3 + j^yin{z 1 z 2 y - ^2 

A\A 2 



1/3n2; 



52.8/ 2 



1.99/ 2 



(A l + A 2 y 

For frontal collisions the fusion radius is still more precisely given by 



RfnsMfm) = (A{ /3 + A^) 



1.908 - 0.0857/n(ZiZ 2 ) + 



3.94 
Z\Z 2 



(19) 



5 Excitation energy 



Assuming a full equilibrated ground state the experimental excitation 
energy E* is the sum of the beam energy E cm in the mass centre frame 
and the Q f US ion value. For the heaviest already investigated experiments, E* xp 
is compared, in table 2, with the energy E£ ar needed to pass classically the 
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outer fusion barrier peak given by the GLDM. The difference between these 
two values is given in the last column. The negative sign corresponds in our 
approach to a tunneling process in the entrance channel. The range of some 
MeV below or above the barrier and the barrier profiles are consistent with 
the very low fusion cross sections and the ability to reach the compound 
nucleus. Let us recall that the fusion barrier height derived from the GLDM 
is systematically higher than the Bass barrier height and the fusion radius is 
smaller [13]. Consequently, for most of the presently used incident energies, 
the fusion is above the barrier for the Bass model and a subbarrier fusion for 
the GLDM. 



6 a decay half-lives 



To determine the potential barriers against a emission, the a decay has 
been viewed [19] as a very asymmetric spontaneous fission within the GLDM. 
The difference between the experimental Q a value or the value predicted by 
the Thomas-Fermi model [14] and the value given by the GLDM has also 
been added at the sphere energy with a linear attenuation factor vanishing at 
the rupture point. Within such an unified fission model the decay constant of 
the parent nucleus is simply defined as A = v P. The assault frequency uq has 
been chosen as 10 20 s _1 . The barrier penetrability P has been calculated within 
the general form of the action integral. Most of the decay path corresponds 
to two-body shapes and the reduced mass approximation is sufficient for the 
inertia. The predicted a decay half-lives agree with the experimental data in 
the whole mass range and also for the known heaviest elements [8,19]. Accurate 
formulae depending only on the mass and charge of the parent nucleus and 
Q a have been extracted from the data. They are recalled here respectively for 
the even(Z)-even(N), even-odd, odd-even and odd-odd nuclei. 



logva |7i/ 2 (s) 
log w [t 1/2 (s) 
login T 1/2 (s) 
logw Ti /2 {s) 



-25.31 - 1.1629A 1/6 v / Z + 
-26.65 - 1.0859A 1/6 v / Z + 
-25.68 - 1.1423A 1/6 v / Z + 
-29.48 - 1.113A 1/6 v / Z + 



1.5864Z 

\/Qa 

1.5848Z 
y/Qa 

1.592Z 

y/Qa 
1.6971Z 



(20) 
(21) 
(22) 
(23) 



The global predictions for the a decay half-lives of the heaviest elements 
within the preceding formulas are given in table 3. If the a decay mode is 
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the main decay mode, then for most of them their half-lives will vary from 
microsecond to some days. Nevertheless, some nuclei ( 319 ' 320 1 26 and 317 124, i.e 
Z around 124-126 and N around 193-194) have a very low Q a and consequently 
a very high a decay half-life. About the same results are obtained within the 
Q a value given by the finite-range Droplet macroscopic model coupled to the 
folded- Yukawa single-particle microscopic model [23]. The calculation of the 
half-lives against symmetric and asymmetric fission is another challenge. 



7 Conclusion 

The barriers preventing the formation of superheavy elements and their 
consecutive a decay have been determined in the quasimolecular shape path 
within a Generalized Liquid Drop Model including the proximity effects be- 
tween nucleons in a neck and the mass and charge asymmetry. The dependence 
on the asymmetry in the entrance channel, on the proton magic number of the 
compound system and on the introduction of the Q value has been studied. 
The shell effects have been determined within the Droplet Model formulae. 

For moderately asymmetric reactions double-hump potential barriers 
stand and fast fission of compact shapes in the outer well is the main exit 
channel. Very asymmetric reactions lead to one hump barriers which can be 
passed only with an energy much higher than the ground state energy of the 
superheavy element. Then, only emission of several neutrons or an a parti- 
cle can stabilize the nuclear system and allows to reach a ground state. The 
formation of superheavy elements via almost symmetric reactions is hardly 
likely. 

In this very heavy mass range the Droplet Model formulae lead to fission 
barrier heights of several MeV or even 10 MeV for nuclei close to the proton 
and neutron magic numbers. The main decay mode being the a decay, the 
predicted half-life depends on the selected theoretical Q a value. Within the 
Q a given by the Thomas-Fermi model, partial a decay half-lives are proposed. 
They agree generally with the known experimental data. 
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Figure captions 

Fig. 1. Selected shape sequence to simulate the transition from two-body 
to one-body shapes and the alpha decay path. The nuclei are spherical when 
they are separated. 

Fig. 2. Macroscopic potential barriers for different reactions leading to 
the 270 110 and 302 1 20 nuclei, r is the distance between the mass centres. The 
vertical bar corresponds to the contact point. 

Fig. 3. Same as figure 2 for 302 1 20 but an adjustment is done to reproduce 
the Thomas-Fermi model Q value. 

Fig. 4. Potential barriers for the A8 Ca+ 2M Pu, A8 Ca+ 248 Cm, 50 Ti+ 2A8 Cm 
reactions. The dashed curve corresponds to the macroscopic barrier. The bar- 
rier given by the dashed and double dotted line incorporates a linear correction 
from the contact point till the sphere to reproduce the Q value. The full line, 
dotted curve and dashed and dotted curve include the shell effects given by 
the Droplet Model assuming respectively a proton magic number of 114, 120 
and 126 and an adjustment to reproduce the Q value. 

Fig. 5. Same as figure 4 but for the reactions 6A Ni+ 208 Pb, 70 Zn+ 208 Pb, 
™Ge+ 208 Pb and 82 Se+ 208 Pb. 

Fig. 6. Same as figure 4 but for the reactions 86 Kr+ 208 Pb and 87 Rb+ 208 Pb. 

Fig. 7. Same as figure 4 but for the reactions 88 Sr+ 208 Pb, lw Cd+ 181 W 
and 10A Ru+ 208 Pb. 

Fig. 8. Symmetric fission barriers for 283 112, 292 116, 295 120 and 311 126. 
The full line, dotted curve and dashed and dotted curve include the shell effects 
given by the Droplet Model assuming respectively a proton magic number of 
114, 120 and 126; the neutron magic numbers being 126 and 184. 

Table captions 

Table 1. Characteristics of the macroscopic potential barriers for the 
heaviest nuclear systems. R m in, Rext, H min and H ext are respectively the posi- 
tions and energies relatively to the sphere of the outer minimum and external 
maximum of the barrier. i? 12 is the distance between the mass centres at the 
contact point. H ext j us is the GLDM external fusion barrier height. 

Table 2. Comparison between the experimental excitation energy E* 
and the classical excitation energy E£ ar at the top of the GLDM barrier. In the 
last column the positive or negative sign of the difference E* xp — E£ ar indicates 
a fusion above or below the barrier. 
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Table 3. Predicted logio[Ti/ 2 (s)] for the heaviest elements and nuclear 
systems versus the charge and mass of the parent nucleus and the Q a value 
given by a recent version of the Thomas- Fermi model [14]. 
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